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A UNIVAC 1108 COMPUTER PROGRAM FOR USE WITH 


WORLDWIDE CLOUD-COVER DISTRIBUTION DATA 

By Kirby D. Kyle 
Manned Spacecraft Center 


SUMMARY 


Worldwide cloud-cover distribution data have been obtained for use in mission 
planning for aircraft and spacecraft missions on which earth resources remote sen- 
sing measurements are made or landmarks are used as navigation aids. To use these 
data, it was necessary to develop the computer, program described, in this report. 

This program, for use with the Uhivac 1108 computer, provides a method of (1) stor- 
ing the cloud- cover statistical data permanently and (2) retrieving the data readily. 


INTRODUCTION 


Worldwide cloud-cover distribution data (ref. 1) were collected and supplied to, 
the NASA Manned ^acecraft Center (MSC). Cloud- cover statistical data are impor- 
tant in mission planning for aircraft and spacecraft missions on which earth resources 
remote sensing is conducted or landmarks are used as navigation aids. These data 
consist of "unconditional" cloud-cover probabilities for specific times of day (desig- 
nated in local standard time (l.s.t.)) over standard areas located by single points; 
each point is specified by a latitude and a longitude designation. "Conditional" proba- 
bilities in terms of 24 hours or 200 nautical miles are also given by the computer 
program. However,, the conditional probabilities can be used at times and distances 
other than 24 hours or 200 nautical miles. The conditional probabilities are pre- 
dictions based on the unconditional (or absolute) probabilities given for a specific 
point at a specific time. The conditional probabilities are for either time into the 
future or distances away from the specific point. It is extremely important that it 
be understood that these statistical data should never, be used for real-time weather 
forecasting, because the applicable atmospheric dynamics are not accounted for in - 
the statistical data. 

The Univac 1108 computer program discussed in this report was provided to 
implement use of the data in preparation for missions, such as the following, that 
include remote sensing experiments. 

1. For each Apollo lunar mission, an earth-orbit contingency mission is planned 
in the event an abort of the lunar portion of the mission, without the necessity for an 



immediate return of the crew, should occur. These contingency earth-orbit missions 
would use the cameras on board the spacecraft for earth resources photography and 
therefore are cloud-cover sensitive. 

2. The Skylab mission to be laimched in April 1973 will carry five sensors for 
earth resources remote sensing measurements. These measurements will be made 
in the visible, infrared, and microwave regions of the spectrum. These sensors will 
be used by the flight crews for a total of 140 days during 1973. Since several hundred 
targets may be established, the long-range mission planning, using cloud-cover sta- 
tistical data, is essential to the success of the mission. In the mission planning, these 
data are correlated with other considerations such as orbital ground track, sun eleva- 
tion angle, and crew activities. 

To use the cloud-cover statistical data, it was necessary to search the data an 
indefinite nmuber of times in order to recall a specific block of the data. The most 
efficient method of conducting this data search was found to be the use of the random 
recall capability of the high-speed-drum storage facility. However, data cannot be 
stor^ permanently on the high-speed-drum facility; therefore, it was necessary to 
provide a capability for rapid transfer of the data from permanent storage to the 
high-speed drum. ■ This requirement was satisfied by the preparation of a computer 
program to transfer the cloud- cover statistical data from punched cards to magnetic 
tape on which the data are stored permanently and are readily transferred to the high- 
speed drum. 

The computer program was written in the FORTRAN V language with the addi- 
tion of the MSC computer system subroutines MREAD, MWRITE, RINIT, RREAD, and 
RWRITE for the Univac 1108 computer. Compatibility with other computer systems 
depends on the availability of a high-speed drum and the ability of the computer system 
to read and produce binary magnetic tapes with a capacity of 38 and 90 words per 
"record" (increment of information). 

Three appendixes are included in this report. Appendixes A and B each contain 
(1) illustrations of the flow diagrams and (2) computer printout listings of the com- 
puter program and subroutines discussed in this report. Appendix C contains a de- 
scription of the magnetic tape storage used with the program. An alternative to the 
method of retrieving the data by use of the high-speed drum is also described. 

Copies of the cards for this computer program are available upon request from 
David E. Pitts, mail code TF321, Manned Spacecraft Center, Houston, Texas 77058. 
Requests for copies of the data on magnetic tape will be honored if the requests are 
accompanied by blank magnetic tapes. 


WORLDWIDE CLOUD -COVER DATA DESCRIPTION 


Two decks of cards contain all the input data necessary to use the worldwide 
cloud- cover statistical data (ref. 1). The first deck contains climatological (map) 
region numbers and boundaries. The second deck contains cloud-cover statistical 
data, compiled for each of 12 months, for each of the 29 map regions. 



Data for the first deck are extracted from the map in figure 1 in the following 
manner. Boundaries of each region fall on even-numbered latitudes and longitudes. 
The area between 70° S and 70° N is divided into 70 swaths at odd-numbered latitudes 
which extend from 0° to 360° eastward from the Greenwich meridian. The areas 
above 70° N and below 70° S require another logic, since one region number defines 
the entire area. Scanning eastward from the Greenwich meridian along each swath, 
the number of the region previously encountered and the value of its terminating longi- 
tude (integral numbers between 0 and 360) are recorded and punched on cards. The 
maximum number of terminating longitudes in one swath is 19. Two cards are used 
to catalog one swath even though data for some swaths do not extend into the second 
card. The card setup is illustrated in figure 2. 

The second deck, containing the cloud-cover statistical data, is illustrated in 
figure 3. Figure 4 illustrates the setup for the individual cards containing the statis- 
tical data. As shown, the data on the cards are organized in three matrices, the first 
being the unconditional probabilities for five cloud categories and eight local times. 

The second matrix contains the 24-hour (temporal) statistical data, and the third : 
matrix contains the 200 -nautical- mile (spatial) statistical data. 

Approximately 75 percent of the cloud-cover data were expressed in tenths of 
sky cover; the remainder was expressed in eighths. Cloud-cover category designation 


is as follows: 



Category 

Tenths 

Eighths 

1 

0 

0 

2 

1, 2, 3 

1,2 

3 ■ 

4, 5 • 

' 3, 4 

4 • - 

6, 7, 8, 9 ; 

5,6,7 

'■■ ■ . . -5; ■ 

10 ; 

8 


Note that the intervals are unequal in size and that the categorization in tenths differs 
only slightly from that in eighths. 

DATA TRANSFER PROGRAM 


The Data Transfer Program transfers only the region number and the termina- 
ting longitude of the map region block to magnetic tape. The transfer is made one swatli 
at a time so that the swath number can be used as an index, for data recall, for an ar- 
ray of points of entry to a high-speed drum. Similarly, only the probabilities of the 
cloud-cover probability data block are transferred. This data block is transferred one 
region at a time (data from five cards) month by month from month 1 to mOnth 12. The 
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month number and the region number are used as indexes for a two-dimensional array 
of points of entry to a high-speed drum. The number denoting the cloud-cover cate- 
gory is not transferred and must be generated by the processing program. The re- 
maining numbers that appear as data on the cards are used to arrange the data card 
deck properly, but are unnecessary for the computational processing and are not 
transferred to magnetic tape. The flow diagram and computer printout listing for the 
Data Transfer Program are given in appendix A. 


DATA RECALL SUBROUTINE SET 


In addition to the Data Transfer Program, two subroutines are used. Subroutine 
CLOUD is the data-manipulating program. Subroutine DRUMST, which is called!,by 
subroutine CLOUD, transfers the tabular cloud-cover statistical data from magnetic 
tape to a high-speed drum. Subroutine DRUMST is called only once for each loading 
of the subroutine. The flow diagrams and computer printout listings for subroutines 
CLOUD and DRUMST are given in appendix B. 

To call subroutine DRUMST, subroutine CLOUD sets all parameters. The call- 
ing statement for subroutine CLOUD is' 


CALL CLOUD(MONTH, HOUR, PLAT, PLONG, UPROB, TPROB, 
SPROB, SCPROB, SCALE, ISCALE) 


where MONTH is the month of the year from 1 to 12; HOUR is the time of day 1. s.t. 
in the range -24. 0 to +48. 0 hours at the point specified by PLAT and PLONG; PLAT is 
the latitude in the range -90. 0° to +90. 0° with north latitudes positive and south lati- 
tudes negative; PLONG is the longitude in the range -180. 0° to +180. 0° with longitudes 
east of the Greenwich meridian positive and longitudes west negative; UPROB, TPROB, 
SPROB, and SCPROB are vectors for storing the probabilities (UPROB (five elements), 
unconditional; TPROB (5x5 array of elements), 24-hour time conditional; SPROB 
(5X5 array of elements), 20O-nautical-mile space conditional; SCPROB (5X5 array of 
elements), scaled conditional (time or space)); SCALE is the scale factor for scaling 
the conditional probabilities, given in hours for time scaling and in nautical miles for 
space scaling; and ISCALE is a control variable which indicates either the type of 
scaling to be done or that no scaling is to be done. That is 


1; time scaling 


ISCALE = 


2; space scaling 


I 3; no scaling 


If ISCALE is not assigned a value of either 1, 2, or 3, it is assigned the value 3 when 
subroutine CLOUD is entered, so that the computer will not call for an error termina- 
tion of the run. 
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The data that must be supplied to subroutine CLOUD are supplied by means of 
the variables MONTH, HOUR, PLAT, PLONG, SCALE, and ISCALE in the calling 
statement. Subroutine CLOUD outputs the unconditional, the 24-hour time conditional, 
and the 200-nautical-mile space conditional probabilities. The scaled conditional 
probabilities also are output if scaling for times other than 24 hours or for distances 
other than 200 nautical miles occurs. 


CONCLUDING REMARKS 


The worldwide cloud- cover distribution data discussed in this report were ob- 
tained for use in mission planning for aircraft and spacecraft missions bn which earth 
resources remote sensing measurements are made or landmarks are used . as naviga- 
tion aids. A computer program, for use with the Univac 1108 computer, was developed 
so that the cloud-cover statistical data can be permanently stored and readily retrieved. 

Manned Spacecraft Center 

National Aeronautics and Space Administration 

Houston, Texas, October 5, 1971 ' 

160-75-03-01-72 
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Figure 4. - Individual punched-card setup for cloud-cover statistical data. 



APPENDIX A 


DATA TRANSFER PROGRAM 



DO loop, index i 







Computer Printout 

I* UIHEnSION MAPC38) |MAP2(3e) .1CL0U0(90) .1ICL00I90I 

2». - - RE.» 1nD I J ^ 

3« ISTAT“U 

<« • U 0 1 0 I " 1 f / u 

S* KEALMbiBUO) ISWTHI I (MAPI Ji iJai tZR) 

6* KEADtS.bOll tSWTH2i ( HAP |J ) , j-29 , 38 ) 

H JSTA1 -ISTAT*! 

6» tiO To (8|H,biAt7,/t7,7«'/ ), jSTAT 

?• 5 HKiTt (A,8J3| IKiISi«Th 3» ISWTH'it ( H AP2 ( J) t J" 1 1 38 ) 

lO» TO 99V 

-II* 6 'ViKITt « 6 , 80 n IK.1SWTH3. ISWTMH,-tMAP2(0) t J*lt38| 

. 12* ISTaT"0 ^ ■ 

13« (iO To B 

n* 7 WHITE (A|8o9| ISTaTi |K 

. IS* GO TO 9V9 - - - 

14* 8 00 9 K*l ,38 

17* 9 HAP2{K) ■MAP(k) 

18 * 1SWTH3" ISwTHl 

I 9* ISaThH- IS«Tm2 

20* 

21 * 10 — CALL MWRITe(l,l,HAP2«nt3a. ISTAT) - 

22* 00 9o K1K8Y«1 ,12 

23* 00 9q KYLE*1 129 

2H* READ(5,5Q2) 11 CLOUD ( I 1 1 , 1 1 ■ 1 . 90 ) - ^ 

29* 3*1 OSTAT-ISIAT*! 

2B» <50 TO (38.3‘<,3St3A,37,37,37i37i37 )*JSTAT 

27* - 35 - WRITE(6i 802) ;K0,KK, 1 1 1 CLOD 1 1 I T • 1 | • t *90 ) - 

28* G0T0999 i 

29* 38 HRITE(8,a03) KO t KK • I I I CLOD 1 1 1 1 * I | • I * 90 ) 

30* ISTAT*0 — ■— ; 

31* GO TO 38 

32* 37 WK|TE(6i8U9i ISTaT,KO,KK 

33* - GO TO 999 — - — 

39* 38 00 8a JK«1 ,9;j 

35* 80 IlCLOOlJKl* KLOUDIJKI 

36* KO-K IHUY 

37* KK-KYLE 

38* 90 caul HWRI rt ( 1 , I , nCLOOl 1 » ,90 i ISTATI 

39* --S00 PORHaTIIX, 12, 19(12, I3»» 

90* SOI P0KHAT(U*12,SU2,IJn 

91* 800 F0RMaT(1X,»E0T MaRK SENSED AlTH l > SWTH ■ t,213,» AND HA 

92*— -IP DATA BEING',/, 2819,/, 10I9> ' 

93* 801 format (U, 'parity ERROR I'rItTEU wITh 1 • »,13,' SWTM ■ *,2I3,» 

99* IAND map data BE 1 NG • , / , 28 I 9 , / , 1 0 I 9 ( 

95*. - 502 - format (6*,l8l2,9(/,6X,18l2> » - 

96* 802 format (IX, »eOT MARK SENSED NITH mONTh ■ ',13,» REGION ■ **I3,» 

97* lANO cloud OAtA BE 1 NG • ,5 ( / ,4X , 1 8 I 3 1 ) 

98*--— 803 -- format (U,»PARItY ERROR WRITTEN PITH MONTH ■ *,13,* REG 1 ON -•--*, 4- 
99* 13,' AND CLOUD DaTA BE i NG • ,S ( / ,8 X , 1 8 I 3 I I 

S3* 900 format (1X,'sooo RON UNLESS PARITY ERROR, *| 

51* .-901- format (IX, 'TAPE NOT LONG ENOUGH TO STORE ALL DATA* BREAK INTO 2 

52* ITAPEs and tNY AGa I N ,*,/, 1 X , 'PROGRAMMER KILLS J0B,*I 

S3* 809 format (IX, 'illegal STATUS*, 19,* ON LOOP *,213» 

S9* _ write I6,9oJ)- ^ 

55* GO TO 998 

S8* 999 WRITE ( 6,9ol | 

57*-- GO TO 997 .. - - _ . . - 

SS* 998 CONTINUE 

59* ENOFILE | 

60* 9.92 rewind 1 - - 

At* CALL EXIT 

62* end 


APPENDIX B 


DATA RECALL SUBROUTINE SET 



Subroutine Set Block Flow Diagram 











Subroutine CLOUD Flow Diagram 

First call only 


Error 

True 

message 





Set 

IREGON = 0 


Calculate random 
access address 


RREAD 38 words 
into MAP array 




Set 

JEFF = 0 



CALL DRUMST 

Data transfer 
from tape to 
high-speed drum 

r 


Set 

IREGON = 15 



Set 

IREGON = 24 



B-4 
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DO loop, index k 



DO loop, index j 







Subroutine CLOUD Computer Printout 


I • 
2* 

3» 

H* 

6* 

/• 

b* 

9* 

1U« 
i I* 

13* 

lb* 

16« 

J/* 


22* 

23* 

2M* 

2b* 

*26* 

27* 

28* 

2V* 

3p* 

il* 

d2« 

i3* 

3H* 

^b* 

36* 

37* 

3b« 

39* 

Hu* 

HI* 

H2# 

H3* 

HH* 

Hb* 

H6* 

H/* 

He* 

HV* 

bU* 

&l* 

b2* 

b3* 

bH* 

bb* 

b6* 

b/* 
bfa* 
bV* 
6U« 
61 • 
62* 
63» 
6H* 
6i>» 
66* 
67« 
66* 
6V» 
/lI» 


s 

•) 

707 


I 


2 

3 


7 

1 1 
12 


n 

13 


17 


18 

19 

30 

33 


31 

35 


3M 


SUbXOuTItHt Cloud I month, HOUR , PL 4T , PL ON 6, UP ROB ,TPROBiSPROB,SCPROB,SC 
lALE.ISCALt) 

DIMENSION MAPA0(70|,M0R6A01>2,79),mAPI 38I,ICL0UDI90),UPB0B(5),TPR0 
IBI5,sli3PM0BIS,5l,SCPR0BtS>5) 

DATA JtFF/0/ 

11- IISCALE iLT. I ,0H, ISCALE «GT. 31 GO TO 5 
GO TO H 

«n1TE16,7U7) ISCALE 
I SCALE*3 


IHL(.0.N"0 

FORMAT!///. • ISCALE NOT ALLOREO. SET ISCALE ■ 3*»» 

IF ( JtFF.LT.bOICALL DRUMS T1 MaPA D , HORGAD , MAP , I CLOUD , JEEF ) 


IF (AUSIPLAT) ,GE« 70.0) GO TO II 
ILAT«IPLAT*71. 01/2,0*0. 5 
NAUUR>MAPA0( ILAT I 


CALL MReau'"aoo'’»”a'’ II>.38,IST*T) 
I lONG*Pi.OnG 

If(ILONG^LT^alILONG“ILON6*3AO 
I ) ( I ST aT .GT . II GO TO 998 

,, I,ST*T)9»0.3.2 

n[) » .-2.38,2 

If (MaP! I I .LT. ilONgI gO TO 7 

INEG0N-MAP(1-1 I 

GO TO 12 

CUNT I NUE 

GO TO 12 

IHEGON-IS 

I F ( PL A T .L T .U . 0 I 1RE60N"29 
.VA0DR6M0RGA0IM0NTN, IREGONI 


call HREAQ(NAUOR,lCLOUOn>,90,ISTATI 
XmOOR-MOUR 

IF (HOUR, lT. 0.01 AhOUK-HOUR*2*I. 
1mOuR-<KHO UR * 0.51/3. 0*1.0 
1FIIH0UR.GT.8I I HOUR* IHOUR-8 
IF ( I STAT .GT. II GO TO 998 
I F ( I STA T1998, 13 . IN 
RnITE(6,700l 
OU 17 I • I , 5 
J« ( I - I I • 1 8* I hour 


UPKOBl I |61CLOUO(JI/IQO. 

amite(6.7oii month, hour, Plat .plong ,iregon,i,uprobiii / 

A R 1 T E ( 6 , 7 (j 3 1 
00 I 9 I • I , b 
OU l« F»1 ,b 

( 1 - 1 I • I 8*9* (A- 1 I 
,i» I 1 - 1 I • 1 8* 1 M* ( K- I I 
TPROR I I ,K I • (CLOUD (LI 
TPKOS( I ,F|6TPR08(1 .KI'/IOO. 

SPRoa ( I ,M •IClOUO ( M I 
SPKUS ( I ,K I -SPROB ( I ,F I / 100. 

AriITEl6,7U2l 1 ,ITPr08(I.KFI.KA-(.5I.I,(SPR0B(I,JJI,JJ-1,5) 

GO TO (30,31,321, (SCALE 
□ U 33 lal ,5 

00 33 K»1 ,S 

SCPR08 ( I , K I -TPROb ( I ,K I 

0 1 /•2M • 0 

IF (SCAuE .GT. 36. nl GO TO 92 
ImEAD« 6H TIHL 
|nLATl*6H MRS* 

AM I TE ( * , 7UM I 1 HE AO ,SC ALE , IHEau I 

GO TO 39 

DO 3S I • 1 , 5 

DO 35 I , 5 

StPMOM ( I , A I 'SPROfi ( I ,K > 

D I V >200 . n 
I ME A0*6H SPA^t 
I HE AO I -AH NM, 

Am 1 t'e 1 * , 7U>( 1 1 hE AO .scale , IHEAO I 

DU 36 K*l ,5 
DO 37 J« 1 , S 
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59 


71* 

72» 

75* 

79* 

75* 

_5B 

IF !K .EO. U) SO TO 56 

SCPROalK . JI9 SCa^9S£PR0B<K (UI/DIV 

IF !scfrob!k*ji .gt. uPROaiUi .and. Scale .gt. divi 

GO TO 57 

SCPR0B!K.J)9|.-SCALE9| |.•SCPROBlK•JI l/DIV 

74* 


IF !SCPROB!K.J) .LT. UPROBIJM GO TO 59 

779 

57 

CONTINUE 

7B9 ' 


SO TO 54 

799 

59 

DO 90 M9|.5 

809 

90 

SCPR0B!K.H)9UPR0BIM) 

8|9 

54 

CONTINUE 

829 


*RITF!4.705I !l.lSCPR0B<l.J>.JPl.5l.t-l.5) 

859 


SO TO 52 

899 

92 

•R|TE!8.704) 


SC TC 


SB* 

• 4* 
87* 
aa* 

89* 

90* 
9l« 
92* 
95« 
9<t« 
95* 
94* 
97* 
98* 
99* 
too* 
V lOl* 
102 * 
105* 
109* 
lOS* 
104* 
107* 
108* 


52 return 

709 FORMAT(////i* SCALE0*iA 4.* COHD I T I ONAt PROS AB I L t T 1 ES* i // i * SEAL 
IE factor is* .F 9i2tA8./i * 61VEN PROBABILITT Or OBSERVINO CATE6 
20RY'«/.' CATE&ORY I 2 3 9 ' 5»»/l 

705 F0RHAT(I4f3X,SF7.2) 

704 F0RHaT(//i* PROBABILITIES ARE NOT TINE CONDITIONAL OVER 54 HOURS* 

1 USE unconditional PROBABILITIES PREVIOUSLY OUTPUT. •> 

998 RNITEI4.750) NADDR. MONTH 

RETURN 

700 FORMATIIHI. • UNCONDITIONAL PROBABITITES OF 0BSERVIN6 CLOUD COVERA 

ISC FOR EACH CLOUD CATCSORY*. //•2Xi'H0NTH HOUR LAT LONS 

2 REStON CATESORY PROBABILITY*./) 

701 F0RMAT(2Xil5i5F9.2,I7tlB.8X.F5.2> 

702 FORmiTII A .*/« . 14. • /• .F4.2 >9F7 . 2 t • /*.I4>* / • tF4 . 2 , 9F7 . 2 • • 

1 /• ) 

705 FORMAT!////. • / • . t 2X . • 29-HOuR T|HE COND I T I ONAL ' . I I X . • / * • I OX . • 200 

INM spatial conditional*. /.* /•.97X,*/*, 97*. •/*,/,* / • . 97X . • / • , 97X . 
2'/'./. * / GIVEN / probability OF OBSERVI 

5ns Category / given / probability of observing catesory /*./i 
9* / category / I 2 5 95/ CATEGORY / 

6l 2 5 95 /*,/.* /*,I0X.*/*,54X.*/*.|6 x.*/*i 

S54X,*/*) 

750 FORMAT! IX. ‘ILLEGAL AOORESS • . 1 1 5 , ♦ MONTH IS *.13,* HOUR IS •»F7.5) 
END 














B-13 





























B-14 













Subroutine DRUMST Computer Printout 


1 • 



' 

H* 


b* 


6* 

lOU 

7* 


b* 

98 

9« 


1 U* 


1 1 • 


12» 

97 

13* 


IH* 


lb* 


16* 

HbO 

17* 


18« 

99 

19* 


d-U* 


21* 


22« 


23* 

H 

29* 


2b* 

b 

26* 


2/. 

7 

2d* 


29* 

a 

3U« 

1 9 

3l» 

lb 

32* 


33* 

18 

JH* 


3b» 


38* 

1 7 

37» 


38* 

1 8 

J9* 

10 

**u* 


91. , 


•«2* 

1 1 

93. 




Hb* 


96* 

29 

*! /• 


Hb» 

2b 

*»9* 


bo* 

27 

bl. 


b2» 

2b 

bj. 


bR. 

3b 

bb* 


bb* 


b/* 

39 

bu» 


b9* 

36 

6u* 


6 1 • 


6^* 

37 

63* 


6h* 

3B 


30 

66* 


6/* 



Subroutine orumst (map*o tHOR6*o,MAP« iclouo , jeff ) 

0 1 mens ION M*PA0(70> .HORS AO I I2.291,HAP(.18l,HAPl08l, |CLOUO< 90» , I ICL 

10 019 01 ^ L__ 

NtE0*33980 

KTLE»0 

Call ri n I T 1 1 BAnoR iNaos I 

00 9fl K»l . 70 

MAPAO ( F 1 ■ I HAOOR* (F- 1 t *38 
FTLE»kTLE*1 

00 97 !■! , 12 

DO 97 J"1 .29 

MoHS/U)<I.«Jl"IbADOR*70»38*Il“l)»9o»29*Tj-r)»9o ' 

IF (NEEO.LE.nAoSI go to 99 

IF (NEEO.OT. NADS. and. FYLE.LT. 1000) GO TO lOU 
1 TE 18. 8S0 ) 

FORM A T ( 1 X . • H I GH SPEED ORUM TOO SMALL TO RUN Jpb . KILL.') 

60 TO 998 
ISTAT«0 
ISTATl-n 
RLAINU 1 
DO I 1 I*| .70 

CALL hRCAUI1.I,MaPM).38.1STAT.LEN) 

JSTAT1"ISTAT1»1 

(.0 TO ta.R.S.T.T.T.T.T.T) .JSTATl 
ARlTe(6.802)MAODR.(HAPl|J)tJ>t.3B) 

GO TO 998 

ARITE<6.803) ISTATl 
GO TO 998 
JSTaT»1STAT»1 

GO TO ( 18 .8. IS. U . 17 , 17. I 7 . 17. 17 ) , JSTAT 
ARITE(6.8U0I 1 .(MAPI J) .jal .38) 

<iO TO 998 

v.R 1 TE I 8,80 1 ) I , (MAP I J) .3*1 .38) 

1 ST A r»o 
GO TO 18 

»RITE(8,8U3) ISTaT 
GO TO 998 
00 10 J«1 ,38 
MAP 1 I U I sMAP ( J I 
NADDR«MAPAD(I) 

IF (lEN .NE. 38) aR|TE(6.8SI) L En , I , JE FF , ( M A P | I J I . Ja I , 38 ) 

CALL RTRIIE(NA0DR,MAPI(1),38.ISTaT1) 

00 2 1 1 « 1 , 1 2 

00 21 J» 1 , 29 

CALL MREA0(l,l,ICLriU0(l).90,ISTAT,LEN) 

JSTAT1«ISTaT1*| 

GO TO (28, 29, 28,27, 27, 27, 27, 27, 27) .JSTATl 
8>ITE(A,7U2) NADOR,(IICLOO(N),Nai,9n) 

GO TO 998 

MX I TE ( 8,703 I ISTATl 
GO TO 998 
JST A r« I ST AT* I 

GO TO <38, 28, js, 38, 37, 37, 37, 37, 37 ) , JSTAT 
IF (I .NE. 12 .AND. J .NE. 29) GO TO 39 
I STAt=U 
GO 10 38 

"RITE(*,7U0) 1,J,( I CLOUD I N ) ,N* I ,90 ) 

Go TO 997 

8RITE(8,7U1) l,J.(ICL0UD(N),Na|,9ni 

ISTAT«n 

GO TO 38 

amITEIA.yJO) ISTAT 
GO To 998 

00 30 K«l,90 

1 1CL0U(K )»ICLUIiD(K ) 

IF (LE'F .NE. 90) m»ITE(8,8SI) LEN , I , j , ( I Icl.OD ( a ) ,K« 1 , 90 ) 

NAL JR»M0RG4ii ( 1 ,J) 



«8* 

21 

CALL RMRiTEINAODRiIICLODU 1 •90t|STAT| ) 



49* 

997 

RE4|.lilD 1 



70. 


JEFF-100 



7 1 . 


RETiibn 



72» 

998 

RE«INO 1 





Call exit 

/.28IH 

,/,iai*i> 

7<(» 

800 

format! IXi'eOf M*RK «EaO on I MjjP |S», 

7b» 

74* 

801 

.ORh^TI 1X.*P^R|TY g^RRoR R^JEO 0*TA 

1 Ir) 

jS' ,/, 

28j*4,/, 10 

77. 

8o2 

FORMAT! |Xi*ILLEGAL AOORESS • • 1 15 1 * ON LOOP*,!],' HAP 

OATA 

IS' ,/,z8l 

78» 


IS,/, 101*41 



7 9* 

803 

FORMAT! IX, « illegal STATUS*,l3t* HAP OATA !S',/,281<| 

,/,IOI<4) 

80* 

851 

format ! 1 X , 'RQROS REAO ERROR LEN ■',!3,' LOOP ■* , 2 1 3 , / , ' 

DATA IS • 

81* 


1 • 5 ! / , 1 8 1 <4 1 1 



82* 

700 

FORMAT! IX, 'EOF MARK READ i ■',13,' J ■',13,* CD ST 

OATA IS', 5!/, 18 

83* 


1 1*41 1 



8M* 

701 

F0RMAT41X, 'PARITY ERROR REAO I •',13,' J ■',13,* CO 

ST DATA IS',S! 

85* 


1/, 181*4) 1 



84* 

702 

F0RHAT!|X, 'ILLEGAL ADDRESS ', 1 1 S , • LOOP',213,' CD ST 

OATA 

IS',S!/,I 

87* 


181*4) 1 



88* 

70 3 

format ! IX ,' illegal status', I 3i* CD ST DXTA IS',S!/, 

ISIH) ) 


89* 


END 
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APPENDIX C 

DESCRIPTION OF MAGNETIC TAPE STORAGE AND 
AN ALTERNATE METHOD FOR DATA RETRIEVAL 

MAGNETIC TAPE STORAGE 


The magnetic tape produced by the Data Transfer Program is a seven-track, 
800-bit/in, binary tape which contains 418 records. The first 70 records of 38 words 
each describe the map regions. The next 348 records of 90 words each contain the 
tabular cloud-cover statistical data. 

In subroutine CLOUD, the vectors MAPAD and MORGAD are used to store the 
high-speed-drum access addresses. The vector MAPAD stores the map region data 
addresses, and MORGAD stores the cloud-cover statistical data by month and region 
addresses. The notation MAPAD.; i = 1, . . , , 70 refers to the ith swath (fig. 1). The 

notation MORGAD.^; j = 1, . . . , 12; k = 1, . . . , 29 refers to the data for the jth month 

and the kth region. In subroutine CLOUD, ILAT corresponds to i, MONTH corre- 
sponds to j, and IREGON corresponds to k. 

ALTERNATE DATA RETRIEVAL METHOD 


An alternate method to the use of the high-speed drum is a magnetic-tape-search 
routine. The magnetic tape record number is stored in the vectors MAPAD and 
MORGAD. The magnetic-tape-search method consists of beginning at the load point 
of the tape and then skipping the necessary number of records to position the tape at 
the required record of information. After the record has been read, the tape is re- 
wound to the load point. 

The following relationships must be used in order to position the tape and read 
the required record. The consecutive numbers 1 to 70 are stored in the 70-word 
vector MAPAD, and the consecutive numbers 71 to 418 are stored in the 348-word 
vector MORGAD. therefore, MAPAD(ILAT) in subroutine CLOUD refers to one of 
the first 70 records on the tape. Similarly, MORGAD(MONTH, IREGON) refers to one 
of the last 348 records on the tape. Then, specifically, MAPAD (33) refers to the 
33rd record, MORGAD(l,2) refers to the 72nd record, MORGAD(7, 28) refers to the 
272nd record, and MORGAD(ll, 29) refers to the 389th record on the tape. The vari- 
able NADDR in subroutine CLOUD contains the number of the record to be read. The 
call to RREAD is replaced by statements to skip NADDR - 1 records and read the next 
record into one of the vectors MAP or ICLOUD, verifying that the first word of the 
record is read into the first word of the vector (i. e. , MAP(l) or ICLOUD(l)). Finally, 
the call to DRUMST and the subroutine DRUMST are eliminated. 
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Retrieval of data by this magnetic-tape-search method is a slow process, and 
better methods do exist. However, the existing computer program, with few changes, 
can be used in the interim during which a more efficient alternate routine can be 
programed. 
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